Decoherence by a chaotic many-spin bath 
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We numerically investigate decoherence of a two-spin system (central system) by a bath of many 
spins 1/2. By carefully adjusting parameters, the dynamical regime of the bath has been varied 
from quantum chaos to regular, while all other dynamical characteristics have been kept practically 
intact. We explicitly demonstrate that for a many-body quantum bath, the onset of quantum chaos 
leads to significantly faster and stronger decoherence compared to an equivalent non-chaotic bath. 
Moreover, the non-diagonal elements of the system's density matrix decay differently for chaotic 
and non-chaotic baths. Therefore, knowledge of the basic parameters of the bath (strength of the 
system-bath interaction, bath's spectral density of states) is not always sufficient, and much finer 
details of the bath's dynamics can strongly affect the decoherence process. 

PACS numbers: 05.45.Pq, 03.65.Yz, 75.10.Nr, 03.67.-a 



Real physical systems are never isolated. Interaction 
of a quantum system with its environment leads to de- 
coherence: the initial pure state of the system quickly 
decays into an incoherent mixture of several states [lj, |2| . 
Modern experiments provide much information about the 
decoherence dynamics of single (or few) ions |3| , Cooper 
pairs 1^1 , or spins [5j , and require a comprehensive theory 
for adequate understanding. Decoherence is also a ma- 
jor obstacle to building a practical quantum computer, 
which, for a wide class of problems, is exponentially more 
efficient than classical computers ||j. Interaction of a 
quantum computer with the bath leads to a fast genera- 
tion of errors, and an accurate theory is needed to find a 
way of controlling this process. 

Decoherence is a complex quantum many-body phe- 
nomenon, and its detailed description is a challenging 
problem. Many theoretical approaches eliminate the 
environment from consideration, approximating its in- 
fluence by suitably chosen operators (deterministic or 
stochastic), and retaining only basic information: the 
strength of the system-bath interaction, characteristic 
energies/times of the bath, etc. |2( Such methods often 
work well, but many situations require detailed account 
of the bath's internal dynamics. Recently, the role of 
quantum chaos 7] in the decoherence process has be- 
come a subject of debate |g, IS llflj - Qualitative semi- 
classical arguments indicate that the chaotic bath (i.e., 
the bath having only a few trivial integrals of motion) is 
"a stronger decoherer" [8| than an integrable bath (i.e., 
the bath possessing a complete set of the integrals of mo- 
tion). Perturbative arguments lead to the opposite an- 
swer |2j. However, these and related works [lfj, eliminate 
the central system from discussion, considering instead a 
static perturbation acting on a bath, and/or treat the 
bath semiclassically, as a particle (or a single large spin) 
with an integrable or a chaotic Hamiltonian. 

Although many valuable insights have been obtained 
in previous work, an important question remains unan- 
swered: is the onset of quantum chaos important for 
the real-world situation when both the system and the 



bath are fundamentally quantum many-body objects 
with non-trivial dynamics? In this paper, we give an affir- 
mative answer to this question. In contrast with previous 
work, we do not replace the system or the bath by per- 
turbation. We go beyond the semiclassical one-body de- 
scription, realistically considering the spin environment 
as many interacting spins 1/2, which have no well-defined 
semiclassical limit. We show that the chaotic bath de- 
coheres the central system stronger and faster than an 
equivalent non-chaotic one, and changes the dynamics of 
the decay of non-diagonal elements of the system's den- 
sity matrix. 

The bath of spins 1/2 (nuclear or electron spins, mag- 
netic impurities) constitutes a major source of decoher- 
ence for nuclear magnetic resonance (NMR) experiments, 
decoherence of pho sphorus spins in Si Jill , spins in mag- 
netic molecules |l2( and quantum dots lljj . Two-level de- 
fects, governing decoherence in Josephson junctions |14J . 
can also be modeled as spins 1/2. Even small coupling 
between the bath spins can make the bath chaotic, and 
we need to understand, at least qualitatively, how this af- 
fects the decoherence process. The dynamics of a system 
decohered by the spin bath is affected by many factors. 
In order to conclusively separate the impact of chaos in 
the bath, and to provide the knowledge needed for more 
complex studies, we need a simple, well-characterized, 
but realistic model. Here, we consider a central system 
of two exchange-coupled spins 1/2, Si and S2, where Si 
interacts with a bath of spins I& (Ik = 1/2, k = 1, . . . N). 
The corresponding Hamiltonian is 



H = JS 1 S 2 + S 1 J2A k I k +H E 



(1) 



where A k are the system-bath coupling constants, and 
Hb is the Hamiltonian of the bath. Similar models 
describe cross-relaxation and double resonance in NMR 
15| , and destruction of the Kondo effect by decoherence 
lfij . Detailed theoretical assessment of specific experi- 
ments requires separate consideration, beyond the scope 
of this paper, but this simplified model captures essential 
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FIG. 1: Time evolution of the elements of the system's 
density matrix p for a chaotic V — 0.04, h = 0.014 (c) 
and a regular V — 0.008, h = 0.014 (r) environment, (a) 
The dynamics of the elements (|| \p\ [[) , (J,f |p| J,f), and 
Im(Tl |p| IT) for b = 0.137. (b) The evolution of the element 
P12 — Re(|| |p| |t) From top to bottom the coupling between 
central system and environment is b — 1.37, 0.137, 0.0683. 
The coupling energy between the two central spins is J = 0.1. 
Everywhere below, the energy and time quantities are dimen- 
sionless. 



physical details of decoherence by the chaotic spin bath. 
Due to similar reasons, we describe the bath by the "spin 
glass shard" model [ijj with the Hamiltonian 



h b = J2 Wf + £ Wl + E h k T k 



k,i 



k 



(2) 



with random T^i and h k ' , uniformly distributed in the 
intervals [—To, Tq] and [0, ho] respectively. This model 
describes the regular-to-chaotic transition in a simple 
and clear way, and permits straightforward control of the 
bath's dynamics |l7j . For small To, the bath is integrable, 
and becomes chaotic for r > T cr ~ h n /(zN) where N is 
the number of bath spins and z is the number of neigh- 
bors coupled via the term Ykil^If ■ Therefore, in real 
baths with large N, the chaotic regime can be relevant 
even for very small couplings Tki- 

We study decoherence by numerically solving the time- 
dependent Schrodinger equation for the wave function 
\^(t)) of the compound system (the central system plus 
bath), using the Hamiltonians Q— (J2J, and considering 
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FIG. 2: The spectral densities of states N(e) vs. energy e for 
regular bath (■) with To = 0.008, ho — 0.014, and for chaotic 
bath (o) with r = 0.04, ho — 0.014. Insets (a) and (b) show 
the level spacing distributions P(s) for chaotic and regular 
baths, respectively. The thick lines show the Wigner-Dyson 
and Poisson distributions for chaotic and non-chaotic baths, 
respectively. 



up to N — 16 bath spins (the results do not change much 
already for N > 10). We use Chebyshev's polynomial 
expansion, in order to work with large Hilbert spaces 
and to study the system's dynamics at extremely long 
times [lj| [l9| . The initial state of the compound system 
is |^(0)) = |</>)|x); where the state of the central system is 
maximally entangled, \(j>) = -^=(| Tl)-| IT)) (singlet), and 
the state of the environment \x) is a superposition of all 
basis states with random coefficients (which corresponds 
e.g. to the bath of nuclear spins at temperatures above 
few tenths of Kelvin) . A wide range of the parameters J, 
ho, To, N, and different sets of Ak have been explored, 
and typical results are presented below. 

It is convenient to describe the system's evolution by 
the reduced density matrix p(t) = Trs|4 f (t))(^ r (t)| where 
Tvb means trace over the bath states. Dynamics of some 
elements of pit) is shown in Fig. ^. Two stages are 
clearly seen: first, the bath rapidly decoheres the sys- 
tem, excites the triplet states, and the system oscillates 
between the singlet and triplet states. Much later, ther- 
malization takes place at much slower rate (note the log 
scale of the time axis). Fig.^i shows the evolution of the 
real part of the non-diagonal element p\% = Re(|| |p| IT) 
for different system-bath couplings, with the coupling 

parametrized by the quantity b — I Y] k _ i A\ ) . For 
every b two curves are shown, corresponding to To > r cr 
(i.e., chaotic bath) and r < T cr (regular bath); for the 
bath here, T cr ~ 0.013. The chaotic bath changes the 
system's evolution both at long time (clearly seen on 
Fig. H 3 ), and a t short times (see below). We verified 
the onset of chaos by calculating the level spacing statis- 
tics P(s) |3, which agrees with the Wigner-Dyson dis- 
tribution for chaotic bath and with Poisson for a regular 
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FIG. 3: The value p^ vs - ^/^ f° r the chaotic environment 
To = 0.04, ho — 0.014 (•), and for the regular environment 
To = 0.008, h = 0,014 (□). The inset shows pf 2 as a func- 
tion of J/b for ho = 1/V2, T = 0.008 (regular bath), for dif- 
ferent sizes TV of the environment and different values 6 of the 
coupling between the central system and the environment: 
TV = 8,6 = 0.518(D), AT = 10,6 = 0.562(c), TV = 12,6 = 
0.683(A), TV = 12,6 = 0.608(H), TV = 12, 6 = 0.965(»), 
TV = 12,6 = 1.365(A). 



bath (insets (a),(b) in Fig- El - It is important that other 
parameters of the bath remain practically intact: the 
large-scale structure of the bath's spectrum (see Eq. |5J) 
is governed by the local fields h^ z , since r cr <C ho. E.g., 
Fig. |21 shows that the spectral density of states is practi- 
cally the same for the regular and the chaotic bath. 

Decoherence can be quantified by the system's entropy, 
concurrence, etc., but particular choice does not affect 
the conclusions. The element p i2 = Re(TJ \p\ IT) is 
particularly suitable for our model: it has an obvious 
physical meaning, and its evolution can be understood 
from the Hamiltonian QJ. The coupling JS1S2 inside 
the central system preserves the initial singlet correla- 
tion between Si and S2, thus steering the value of pi 2 
towards —1/2. The system-bath coupling Si ^A^lk en- 
tangles the spin Si with the bath and destroys the cor- 
relations between Si and S2, thus leading p\% towards 
zero. Competition between the two tendencies deter- 
mines the value of P12 at t — > 00. Inset in Fig. |3| shows 
that P12 = Pi2(t — > 00) is determined by the single ratio 
J/b (where b 2 = ^-A^), independently of the size of the 
bath TV and particular values of Ak . But the bath's inter- 
nal dynamics noticeably affects the dependence p\ 2 {J /b) : 
as Fig.^Jj shows. Fig.|3]presents the results of many sim- 
ulations, comparing the curves Pi 2 (.j r /b) for chaotic and 
regular environments. The chaotic bath, for the same 
value of J/b, is more efficient in steering p\ 2 towards zero, 
i.e. in breaking the correlations between Si and S 2 . 

But the most obvious difference between the chaotic 
and the regular baths emerges at short times, t < 100- 
300 (Fig. HJ. For J>i (small system-bath coupling), 
Pi2(t) oscillates with the frequency us ~ J, mirroring the 



-0.45 



-0.50 




-0.490 



-0.495 



-0.500 




100 100 200 300 

time 



FIG. 4: Short-time evolution of the off-diagonal element 
pi2(i) in the case of a chaotic environment To = 0.04, 
ho = 0.014 (solid curve) and in the case of a non-chaotic 
environment To = 0.008, ho = 0,014 (dashed curve). The 
coupling between the two central spins is J = 0.4. The values 
of the interaction between the central system and its environ- 
ment are 6 = 0.2 (a), and 6 = 0.0683 (b). For 6 = 0.2 (a), the 
oscillations for the chaotic and the regular bath are almost 
identical, so that the solid and the dashed curves coincide. 



quantum oscillations of the central system between the 
singlet and triplet states. For b larger than the spectral 
width W of the bath (here, W ~ 0.1, see Fig. E), the 
oscillations of px%(t) are identical for the chaotic and the 
regular baths (Fig.QJi). The envelope of the oscillations is 
Gaussian, i.e. p^^) = a + /3exp (— i 2 /T s 2 ) where a and 
/? are constants, and T s is the decay time. However, when 
b becomes smaller than W , prominent differences appear 
(Fig. ^Js) . For the regular bath, the decay remains Gaus- 
sian, but the chaotic bath leads to the exponential decay, 
with the envelope p c ^{t) = a! +/3' exp (-t/T s ). The de- 
cay time T s also becomes different, see Fig. |5l where T s 
is plotted as a function of 1/6. The values of T s were 
determined from the least-square fits of Pi 2 v (i) to both 
Gaussian and exponential forms; both forms give similar 
T s . For b < W (large 1/6), the chaotic-bath decoherence 
is faster by a factor of 2-2.5 [20|. 

It is important to note that the curves pi 2 (i) for the 
regular bath are insensitive to To, but the drastic differ- 
ence emerges as soon as Tq exceeds r cr , when the bath 
becomes chaotic. 

In our model, the central system can not be consid- 
ered as a Hamiltonian perturbation acting on the bath: 
Si can not be replaced by a fictituous magnetic field 
since the intra-system and the system-bath couplings 
are isotropic |2lJ. Also, our bath has no semiclassi- 
cal limit. Nonetheless, there is a striking analogy be- 
tween our results and the Loschmidt echo decay in semi- 
classical systems [lfj. Following Ref. |8|, if the central 
system could be replaced by a perturbation A of the 
bath's internal Hamiltonian J2J, then the states | fl) 
and I jt) of the central system would correspond to dif- 
ferent perturbations A' and A", which would produce 
different bath states \x' t ) = exp (— i[Hb + A']t)\xo) and 
\x") = exp(-i[H B + A"}t)\ X o)- The strength of the 
system-bath interaction 6 then would correspond to the 
magnitude of ||A||, and the matrix element (fj. \p\ 
would correspond to the overlap (x't\Xt)- R i s known 
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FIG. 5: Decoherence time T 3 as a function of 1/6. Circles 
denote the case of a chaotic environment (Fo = 0.04, ho — 
0.014), squares denote the case of a non-chaotic environment 
(F = 0.008, ho = 0,014). Fo obtain T s we fitted ptT{t) to 
Gaussian (empty symbols) and exponential (solid symbols) 
forms. Inset shows the same data for 1/6 > 10 in a semi-log 
scale. 



that the quantity F(t) = \{x't\Xt)\ 2 (called Loschmidt 
echo) exhibits Gaussian decay when ||A|| is larger than 
the bath's spectral width W, and our results for /312(f) at 
b > W also show Gaussian decay. At ||A|| < W (in Lya- 
punov's regime), the Loschmidt echo of the chaotic bath 
decays exponentially with the rate independent of ||A||, 
while for the regular bath the decay is Gaussian. Our 
simulations give the same picture, with the decay time 
T s almost independent of b for chaotic bath (T s changes 
by only - 20% for 0.005 < b < 0.1). So, it is likely 
that the b < W regime of decoherence corresponds to 
the Lyapunov's regime of the bath, in s pite of the fact 
that our bath has no semiclassical analog |22j . The study 
of the Lyapunov's exponents is an interesting problem for 
further research. 

Summarizing, wc compare decoherence of a two-spin 
system by regular and chaotic spin baths. We go beyond 
the standard one-body semiclassical description, consid- 
ering environments of many spins 1/2. We do not replace 
the system by a perturbation acting on the bath, thus go- 
ing beyond the Loschmidt echo studies. At t — > 00, the 
chaotic bath leads to smaller values of the system's den- 
sity matrix element (fl \p\ J.T) than the regular bath, 
i.e. at long times the chaotic bath decoheres the system 
more efficiently. At short times, the chaotic bath leads to 
faster decay of quantum oscillations in the system, and 
changes the form of the decay from Gaussian to exponen- 
tial. Therefore, the onset of chaos in the bath drastically 
changes the decoherence dynamics. Also, based on the 
analogy with the Loschmidt echo studies, we give argu- 
ments that the chaotic bath is in the Lyapunov's regime. 
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